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ABSTRACT 

At radio wavelengths, images of the compact radio source Sagittarius A* 
(Sgr A*) in the Galactic Center are scatter broadened with a dependence 
due to an intervening ionized medium. We present VLBI observations of 
Sgr A* at 86 GHz using a six station array including the VLBA antennas at 
Pie Town, Fort Davis and Los Alamos, the 12 m antenna at Kitt Peak and 
the millimeter arrays at Hat Creek and Owens Valley. To avoid systematic 
errors due to imperfect antenna calibration, the data were modeled using 
interferometric closure information. The data are best modeled by a circular 
Gaussian brightness distribution of FWHM 0.18 ± 0.02 mas. The data are also 
shown to be consistent with an elliptical model corresponding to the scattering 
of a point source. The source structure in the N-S direction, which is less well 
determined than in the E-W direction due to the limited N-S {u, v) coverage 
of the array, is constrained to be less than 0.27 mas by these measurements. 
These results are consistent with extrapolations of intrinsic structure estimates 
obtained with VLBI at 7 mm wavelength assuming the intrinsic size of Sgr A* 
has a greater dependence than A*^'^ with wavelength. 

Subject headings: Galaxy: center — galaxies: individual(Sagittarius A*) — 
techniques: interferometric — scattering 
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1. Introduction 

Sgr A* is the compact radio source in the center of the Galaxy and is beheved to be 
powered, as are extra-Galactic AGN, via accretion onto a massive black hole. Evidence for 
this scenario solidified recently with the stellar proper motion studies of Eckart & Genzel 
(1997) and Ghez et al. (1998) which show velocity dispersions of ~ 1400 km s"^ near the 
Galactic center implying a 2.6 x IO^Mq central mass within a 10~^ pc^ volume. Upper limits 
on the proper motion of Sgr A* itself (Backer & Sramek 1999, Reid et al. 1999), constrain 
its mass to be greater than ~ 2 x 10"^ M©. These results would seem to rule out stellar type 
phenomena as models for the radiative output from Sgr A* which spans radio to x-ray 
wavelengths. Estimates of the luminosity of Sgr A* (Zylka et al. 1992) though, are only a 
small fraction (~ 10^^) of the Eddington limit for a 2.6 x lO^M© black hole (Quataert et 
al. 1999). Physical models that reproduce the observed Sgr A* emission spectrum in the 
context of mass estimates above include scaled AGN models (Falcke & Biermann 1999), 
advection dominated spherical accretion flows (ADAF) (Narayan, Yi & Mahadevan 1995) 
and cyclo-synchrotron emission from spherical accretion (Meha 1992). 

Direct measurement of the intrinsic size of Sgr A* is of great importance in the 
discrimination of different emission mechanisms thought to operate in this source. In the 
radio bands, such measurements are complicated by scattering through the inhomogeneous 
ionized ISM which broadens images. Past VLBI observations ranging in wavelength from 
6 cm to 1.3 cm, observe an elliptical source with axial ratio 0.53 oriented E-W whose 
dimensions vary with as expected for scattering of an unresolved source by electron 
density irregularities (Lo et al. 1985, Jauncy et al. 1989, Marcaide et al. 1992, Alberdi et 
al. 1993). Similar elliptical scattering disks are observed for OH masers located within 25 
arc minutes of Sgr A* (Frail et al. 1994). Such non-circular scattering implies an anisotropy 
in the underlying density inhomogeneities that may arise due to ordered magnetic fields in 
the region. 

Using 7 mm wavelength VLBI, Lo et al. (1998) reported a departure from this axial 
ratio of the scattering ellipse, a result consistent with independent measurements by Bower 
& Backer (1998). The apparent departure from the scattering, if confirmed, would be strong 
evidence for intrinsic structure in Sgr A*. A possible elongation of Sgr A* at 7 mm was 
first reported by Krichbaum et al. (1993) who modeled Sgr A* with two VLBI components 
separated by 3 mas along a PA of —25°. Later 7 mm VLBI work by Backer et al. (1993) 
and Bower and Backer (1998) found no evidence for such structure. Lo et al. (1998) 
interpret their 7 mm result as intrinsic structure consisting of an elongation of 0.44 ± 0.09 
mas with position angle of —10°. The reported statistical significance of this result is 5 
standard deviations, but there are additional systematic errors in the calibration (Bower et 
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al. 1999). VLBI observations of Sgr A* at 3.5 mm wavelength, which are much less affected 
by interstellar scattering, can provide important verification of the 7 mm results. Intrinsic 
structure of Sgr A*, if visible at 7 mm, should be readily observable at shorter wavelengths 
as long as its size does not decrease faster than A^. 

Previous size estimates for Sgr A* from VLBI at 3.5 mm wavelength have been 
obtained using arrays of two (Rogers et al. 1994) and three (Krichbaum et al. 1998) 
antennas. The sparse data sets of Rogers et al. (1994) and Krichbaum et al. (1998) were 
modeled with circular Gaussian brightness distributions and yielded sizes of 0.15 ±0.05 mas 
and 0.19 ± 0.03 mas respectively, consistent with a extrapolation of the scattering size at 
this wavelength. There remain some differences in the literature as to the exact power law 
which best characterizes the relation between scattering size and wavelength (Backer et al. 
1993, Rogers et al. 1994, Krichbaum et al. 1998). Differences in the number and weighting 
of past size measurements are largely responsible for the variation in fitted relations. In this 
paper, we will adopt the relations cited by Lo et al. (1998) which yield a size of 1.43 A^ mas 
for the major axis and 0.76 A^ mas for the minor axis where A is measured in cm. These 
result from nearly simultaneous observations reduced in a uniform manner. 

The limited baseline coverage and sensitive dependence on a-priori antenna calibration 
of previous arrays have not yet allowed a definitive estimate of intrinsic structure at 
shorter wavelengths, though a detection of Sgr A* at 1.3 mm wavelength has been reported 
(Krichbaum et al. 1998). This paper reports on new VLBI observations made at 3.5 mm 
wavelength with an array of six antennas that allowed use of closure quantities to model 
the Sgr A* structure. We explore two different methods of model fitting using the closure 
amplitude information. Closure amplitudes and phases provide structural information 
without the need to cahbrate the VLBI array (Readhead et al. 1980). 

2. Observations 

Observations of Sgr A* (a(J2000) = 17''45™40^045, 5(J2000) = -29°00'27.9"; Rogers 
et al. 1994) at 86 GHz were made by the Coordinated Millimeter VLBI ArrayQ (CMVA) 
in April 1999. Fringes were obtained to the KittPeak 12m antenna (K), the Owens Valley 
mm- Array (O), the BIMA array at Hat Creek (H), and the NRAO[| VLBA sites at Los 
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Alamos (L), Pie Town (P), and Fort Davis (F), where the symbols to be used in this paper 
for each site are in parentheses. Sgr A* and continuum cahbrators NRAO 530, OV-236 were 
observed on both 16 and 18 April 1999, but technical problems at the VLBA sites reduced 
the array to just 3 antennas on the first day. The data presented here were observed during 
the second day in which Sgr A* was detected on baselines to all sites. A total bandwidth 
of 56 MHz consisting of 14 separate 4 MHz channels was observed for 360 second duration 
scans. SiO masers towards VX Sgr were observed at both 43 GHz and 86 GHz for pointing 
and calibration respectively. The data were processed on the Mark HIA correlator at 
Haystack Observatory with an integration time of 1 second. 

3. A-Priori Calibration 

Accurate a-priori amplitude calibration of VLBI visibilities relies completely on 
the successful determination of telescope sensitivity as a function of time and elevation 
using previously established gain information. At the low elevations required for Sgr A* 
observations, mm-wave dishes such as the Kitt Peak 12 m can be accurately characterized, 
but gains for the VLBA antennas, which are not yet optimized for 86 GHz observations, 
are difficult to determine. Phasing at the mm arrays also adds a degree of gain variation 
that depends sensitively on time variable atmospheric conditions. Two independent a-priori 
calibration techniques were applied to the array. 

The SEFD (System Equivalent Flux Density) was first determined for each antenna as a 
function of time using system temperature (Tsys) measurements and previously established 
gain-elevation curves. VLBA antennas use internal noise sources to determine Tsys, and 
these values were corrected for atmospheric opacity efi^ects. The mm arrays and Kitt Peak 
use a vane absorber method to measure Tsys and no opacity correction is required. Explicit 
assumptions in this type of array calibration are perfect pointing and stable gain curves at 
all sites. A faulty calibration source prevented Tsys from being recorded at the VLBA Fort 
Davis site. 

A second calibration method made use of total power spectra obtained 
by auto correlating VLBI data on the 86 GHz SiO maser source VX Sgr 
(q;(J2000) = 18'^08"^04^04, 5(J2000) = -22°13'26.90"; Wright et al. 1990) to 
generate 112 channel spectra in a 4MHz bandpass. A template to the maser line structure 
from the Kitt Peak auto correlations was calibrated assuming a fiat gain curve and a 
sensitivity of 0.025 K/Jy. This template spectrum was then fitted to total power spectra 
extracted from all VX Sgr observations at all array sites to determine SEFD of each 
antenna as a function of time. This method of calibration assumes perfect pointing only for 
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the template spectrum, pointing errors on all other VX Sgr observations will be corrected 
relative to the template scan. At BIMA, a phase calibration tone leaking into the bandpass 
prevented calibration using the spectral template technique. 

Fig. m shows the results of these two calibrations for all sites. At four sites (Kitt Peak, 
OVRO, Pie Town and Los Alamos) both the spectral template (filled circles) and gain 
curve (solid line) calibrations were used. The Kitt Peak and OVRO sites show generally 
good agreement between the two methods although it appears that the gain curve at Kitt 
Peak may have more elevation dependence than originally assumed. The gain curve derived 
SEFD at Kitt Peak rises more slowly at lower elevations than the VX Sgr derived values. 
For the two VLBA sites, correspondence between the two calibrations is distinctly inferior. 
In the elevation range 20° to 30°, the gain curve calibration is consistently below the VX Sgr 
SEFD values by as much as a factor of 4 in the case of Los Alamos, and by a factor of 2.5 
for Pie Town. The difference between the calibration schemes is probably due to pointing 
errors at the VLBA sites which are given, along with their effects on VLBI amplitudes, 
in Table |I|. The loss factors given are 1 a values, so the observed scatter in the spectral 
template results is entirely consistent with pointing errors as well. At 86 GHz, the VLBA 
pointing is affected by known systematic pointing offsets at the 20 arc second level due to 
ripple in the azimuth track (Dhawan, private communication). Both a-priori calibration 
methods thus do not compensate for pointing errors on Sgr A* observations at the VLBA 
antennas. In order to avoid these difficulties in calibration, which are likely not tractable at 
3.5 mm, we chose to use only closure amplitudes in our analysis which are not susceptible 
to calibration errors and are determined on observations of Sgr A* itself. 



4. Data Analysis 

The Sgr A* data was analyzed using the following procedure: 

1] Fringes were obtained at all sites on the calibration sources to determine the 
calibration phases which allow coherent bandwidth synthesis of the baseband channels. In 
addition the calibration sources were used to estimate the clock delays and clock rates for 
each station. 

2] A fringe search was conducted on Sgr A* to refine the delays and rates using 
the Haystack Observatory Postprocessing System (HOPS) software. This step used 
interpolations of delay and rate from adjacent calibrator scans as well as strong Sgr A* 
detections on sensitive baselines to constrain searches on weaker baselines. This method 
works even when there are very weak fringes on some baselines, as long as there are 
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sufRciently strong fringes on enough baselines to reduce the search space (e.g. Alef & 
Porcas 1986). The data were also imported into the NRAO AIPS software package where 
independent fringe searches resulted in consistent delay and rate results. 

3] The noise bias corrected incoherent average of 10 second data segments was 
computed for each basehne using the refined clock delays from the constrained fringe search. 
These data segments are short enough to avoid significant signal loss due to atmospheric 
path fluctuations. The incoherent average (a) (Rogers, Doeleman, Moran 1995) is 



where is the complex correlation from the coherent integration of each segment and Gc is 
the noise bias 

a,= - {2BT)--^ (2) 

where B = 56 MHz 

T = coherent integration =10 sec 

L = loss factor ~ 0.55 for digital MK3 1-bit correlation including fringe rotator loss 
M = number of segments 

4] Closure phases, the sum of interferometric phases around a closed loop of three 
basehnes (Jennison 1958, Rogers et al. 1974), were estimated for each 360 second VLBI 
recording scan. The complex valued bispectrum or triple product (Rogers, Doeleman & 
Moran 1995), whose argument is the closure phase, was calculated for each 10 second data 
segment and averaged to a full scan length. This avoided any coherence losses due to 
phase fluctuations over the length of an entire 360 second VLBI observation. Deviations of 
the closure phase from zero on any triangle would indicate an asymmetric component in 
the structure of Sgr A*, zero closure phases are the signature of a brightness distribution 
symmetric about the origin. For a brightness distribution dominated by a central compact 
feature, closure phases are not expected to change significantly on a 360 second time scale. 
Closure phase plots in Bower & Backer (1998) show that even for a 2 component model 
of Sgr A* at 7 mm, closure phases vary by < 5° in a six minute period. In the low signal 
to noise (SNR) regime, the bispcctral SNR is proportional to the product of the SNRs of 
each contributing baseline and it can be shown (Rogers, Doeleman & Moran 1995) that this 
SNR is greater than that of the closure phase alone. We employ a bispcctral SNR cutoff 
which ensures an SNR in each 10 second segment on each baseline of approximately 1.6. 
This corresponds to a 15° standard deviation in the closure phase. 
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The measured closure phases are consistent with zero values. Several representative 
closure phase triangles whose high SNR enabled us to track their values over multiple scans 
are listed here with their measured standard deviations : FKP = 24° ± 12, FKL = 10° ± 23, 
KOL = 11° ± 22, FKH = 3° ± 20, FKO = -6.3° ± 27. Subsequent modeling was restricted 
to symmetric brightness distributions. 

5] A circular or elliptical Gaussian is fit to the amplitudes allowing separate station- 
based gains for each 360 second scan. This is similar to the commonly used self-calibration 
technique for imaging interferometric data except that self-calibration cannot search 
all possible brightness distributions. Instead, self-calibration converges iteratively on a 
consistent set of complex antenna gains and an image by using a deconvolution method 
such as CLEAN. The resulting image is not unique (Schwarz 1978) due to the iterative 
nature of the imaging process. Alternatively, the method described below performs an 
exhaustive search over all models in a restricted parameter space and guarantees that the 
result is the most likely. 

The best fit model is found by minimizing x^- For a weighted least squares, 

- E^^.W - m,,eS^We^^W}' (3) 

where my(i)is the model at time t for baseline ij 

Qij (t) is the observed incoherent average of data segments 
g9i(t) is gaij^ for station i for each scan at time t 

and Wijit) = l/aij^{t) 

where cry^(^) is the estimated variance in ay(t) 

The standard deviation in each noise bias corrected amplitude estimate in correlation 
units is 

a^aM-'^{{s' + l)/s')'' (4) 

where s is the coherent SNR for each segment. Monte Carlo simulations show this 
approximation to hold for a coherent SNR =a/af. > 0.5. For example, a correlation of 
3 X 10~^ has a SNR = 0.5; with 10s segments and M = 36, a ~ 1.5 x 10~^. Baselines with 
correlations below 3 x 10~^ were not used in the model fit. 

The summation is carried over all available baselines and scans. For each scan and 
each model the station gains are estimated by the weighted least squares (Bevington & 
Robinson 1992, see also Scharf 1991 for matrix notation) solution of 



y — Ax + n 



(5) 
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given by 



X = 




(6) 



where 



X 



[go, ■■.,gL-i]'^ = gain vector for L sites which participate in scan 



w = 


diagonal weight matrix 


y = 


data vector 


A = 


steering (or curvature) matrix 


n = 


noise or error vector 


T = 


the transpose operator 


X = 


weighted least squares estimate of x 


-1 _ 


-- the inverse operator 



The nth element of the data vector is yn = ln(a„/m„) where n is the baseline number. 
The design matrix elements are 

Ank = 1 when k = i or k = j 
= otherwise 

where k is the station number. The use of logarithms linearizes the estimation of the 
gains. The search for a minimum using Eq. ^ requires separate gain solutions for each 
model tested since the gain solutions depend on the model. A 1-D Monte Carlo search is 
made for the circular Gaussian with minimum x^, while an elliptical Gaussian requires a 
3-D search through parameter space. Pearson (1995) discusses the least squares fitting of 
Gaussian models to VLBI data and the estimation of errors. Under the assumption of a 
Gaussian distribution of errors, the weighted least squares approach, which simultaneously 
solves for the model and the station gains, is equivalent to the maximum likelihood solution. 
Using this method we obtained the best fit circular and elliptical Gaussian models shown 
in Table ^ (models A and B). 

The resulting least squared station gains are directly compared with the a-priori 
results in Fig. 0, represented by crosses. In most cases, the modeled station SEFD values 
are less than spectral template values at comparable elevations. This may be due to the 
proportionality of the modeled SEFDs on the total flux density of Sgr A* which we have 
assumed to be 1.4 Jy based on BIMA observations taken prior to the VLBI observations 
(Wright, private communication). If the flux of Sgr A* during our observing epoch was 
higher, the derived station gains from the modeling would increase. 

Figure ^ compares the gain corrected visibility amplitudes with a model for the best fit 
circular Gaussian. The same data is recast as a plot of correlated flux density as a function 
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of baseline length in Figure ^. The fit to an elliptical Gaussian is not well determined in the 
N-S direction owing to the reduced {u,v) coverage in this direction (see Figure |^). Errors 
in the circular and elliptical Gaussian models were estimated from the size of the region in 
parameter space corresponding to Ax^ < 1. In Figure |^ we have plotted & function 

of the FWHM of the circular Gaussian model. We explored a range of coherent averaging 
times from 4 to 60 seconds and found very little change in our results from 4 to 20 seconds. 
At 60 seconds Xu"^ w^is over 1.5 and the estimated error in the size of the best fit circular 
Gaussian more than doubled. 

We also modeled the data using closure amplitude quadrilaterals which are insensitive 
to station gains. The closure amplitude is defined as : 

^ «ik «jl 

a ratio of complex amplitudes around a closed quadrilateral loop of baselines. The number 
of independent closure amplitudes is equal to the number of baselines between all sites minus 
the number of sites so that for 6 sites there are 9 independent closure amplitudes. Cornwell 
(1995) discusses the equivalence of self-calibration and modeling of closure quantities in 
determining source structure. However, using from fitting the closure quadrilaterals to 
estimate the errors in source parameters is not straightforward as the quadrilaterals are 
not all independent. Furthermore, in the low SNR case, the formal error for each closure 
amplitude is difficult to determine (see e.g. Trotter, Moran, Rodriguez 1998, eqns. 4 and 5). 

Based upon VLBA measurements between 6 cm and 7 mm wavelength, Lo et al. (1998) 
have suggested that the 7 mm results show deviation of the minor axis size from the 
dependence implying an intrinsic size of 0.44±0.1 mas along PA=— 10°. Since the scattering 
angles are smaller at 3.5 mm wavelength, it is important to compare the constraints set by 
the 3.5 mm data to reasonable extrapolations from the Lo et al. (1998) inferred intrinsic 
structure. Model C constitutes pure interstellar scattering of a point source. Model D 
is the extrapolated apparent source size based on an intrinsic source size that varies as A°'^ 
(Lo et al. 1999). Model E represents the case in which the 7 mm derived intrinsic structure 
does not vary with wavelength. Estimates of apparent source size at 3.5 mm and 1.3 mm 
(Rogers et al. 1994, Krichbaum et al. 1998) are consistent with some sort of wavelength 
dependent intrinsic structure from 7 mm to 1.3 mm, but these results depend sensitively on 
a-priori calibration. The exact behavior of the Sgr A* intrinsic structure as a function of 
frequency is thus largely uncertain with estimates for a power law dependence (A°) ranging 
from a = 0.7 to 1.9. Table ^ contains characteristics and fit results for models C, D and E. 
For all elliptical models, the reported PA reflects the major axis orientation. 
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5. Discussion 

The 3.5 mm observations are consistent within 2 cr of the elliptical scattering model 
of a point source. They are also consistent at the 3 a level with a model combining the 
elliptical scattering disk with an intrinsic structure extrapolated from the Lo et al. (1998) 
7 mm estimates assuming a A°'^ dependence. Increasing the spectral index of this power 
law, and consequently decreasing the intrinsic size, improves agreement with the 3.5 mm 
data as the observed size approaches the scattering size. To further illustrate the level of 
agreement of the 3.5 mm data with all models, we show in Figure ^ a closure amplitude for 
the Fort Davis, Kitt Peak, OVRO, Hat Creek quadrangle (C*pkOh) ^long with the same 
closure amplitude values expected for models A through E. Model E, shown as the sohd 
line in the lower panel of Figure ^ provides the poorest fit to the data. The dotted hues 
show the effects on model E if the N-S intrinsic structure from Lo et al. (1998) is reduced 
or increased by their cited 1 a errors. Other quadrilaterals also show similar inconsistencies 
with model E. 

It has been customary in VLBI studies of Sgr A* to model the brightness distribution 
as an elliptical Gaussian given the scattering morphology seen at longer wavelengths. At 
shorter wavelengths (A < 7 mm), however, it is not clear that the data warrant a model 
more complex than a circular Gaussian. Model A in Table ^ has the lowest of any model 
tested, but it also contains two extra model parameters relative to Model B. The F-^-test 
(Bevington & Robinson 1992) is the proper method to determine the statistical significance 
of a drop in when additional model parameters are added. The result of such an analysis 
is based on trying to fit both circular and elliptical Gaussian models to pure noise and 
calculating the difference in x^- Given the number of data points in our observations, we 
find that a Ax^ = 5 between the two models occurs at least 12% of the time when fitting 
pure noise. The better fit of model A, therefore, is only significant at the 1.5ct level, and we 
conclude that the data presented here are adequately fit by a circular Gaussian. A similar 
test has not yet been performed for previous 7 mm wavelength data sets (Lo et al. 1998, 
Bower & Backer 1998), and such an analysis would be useful in determining the significance 
of their results. 

Looking for departures from extrapolated scattering sizes due to intrinsic structure is 
difficult with current arrays. Based on contours of Ax^ = 14 we can limit the elongation in 
a position angle of —10° to less than 0.27 mas. At 7 mm the effects of intrinsic structure 
are small and imaging can be influenced by calibration of the VLBA antennas. Previously 
reported measurements at 3.5 mm by Rogers et al. (1994) and at 1.4 mm by and Krichbaum 
et al. (1998) have insufficient baselines to be able to rely on closure and are consequently 
limited by calibration errors. Closure phases derived by summing interferometric phase over 
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station triplets are generally consistent with zero value. This constrains possible brightness 
distributions to be symmetric about the origin, a characteristic shared by all the models 
tested here. Adding secondary components to the Sgr A* structure produce marginally 
better fits to the closure phase data, but not enough to justify the additional degrees of 
freedom. The closure amplitude data presented here show no evidence of asymmetry. 

Assuming the scattering and intrinsic structure add in quadrature, limits on the 
intrinsic size for Sgr A* in position angles of —10° and 80° are 0.25 mas and 0.13 mas 
respectively. For a 1.4 Jy source, these limits correspond to a minimum brightness 
temperature of 7 x lO^K, nearly half the value derived in Rogers et al. (1994). The essential 
difference between the two results is the assumed structure, which Rogers et al. (1994) 
models as a circular Gaussian. The limits derived here assume an elliptical model allowing 
for an extension in the N-S direction and consequent increase in size. 

Ozel et al. (2000) have recently computed new ADAF models incorporating generalized 
electron energy distributions that reproduce the observed Sgr A* emission spectrum. Their 
size estimate at 7 mm closely matches the size of the N-S extension seen by Lo et al. (1998), 
but is much larger than the upper limit to the intrinsic size set in the E-W direction. The 
spherical ADAF models may have limited applicability to Sgr A* if an intrinsic ellipticity 
or asymmetry is verified at shorter wavelengths. At 3.5 mm, the ADAF model (Ozel et al. 
2000) predicts an intrinsic size of 0.29 mas, well over the 3 a upper limit in the N-S direction 
derived here (Figure |^). In the AGN model of Falcke & Biermann (1999), the estimate of 
intrinsic size for Sgr A* at 7 mm is used to constrain parameters of the relativistic jet, and 
the emission region size varies as A°'^^^° where Xo is a factor related to the jet orientation. 
An extrapolation to 3.5 mm gives a size of 0.24 mas, larger than the 3.5 mm E-W size limit 
but consistent with N-S limits. 



6. Summary 

We have used amplitude closure information from a 86 GHz VLBI array to determine 
a size estimate of Sgr A* without the need to calibrate the array with a-priori information. 
With an assumption of elliptical Gaussian structure, the data are consistent with interstellar 
scattering of an unresolved source. Our results also do not exclude extrapolations from 
estimates of intrinsic size at 7 mm wavelength as long as the intrinsic size evolves faster 
than A°'^. A best fit elliptical model does exhibit an elongation in the N-S direction, 
but the improvement in fit over a circular Gaussian model was found to be of marginal 
statistical significance. The N-S {u, v) coverage of the VLBI array is sufficient to constrain 
an intrinsic extension in the —10° position angle to be less than 0.25 mas. Emission models 
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that reproduce the spectrum of Sgr A* predict sizes for Sgr A* that are at or above this 
size hmit. 
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Table 1: Signal Loss Due to Pointing Errors at VLB A Sites 



Antenna 


rms pointing 


(arcsec) loss factor 1 


loss factor 2 


loss factor 3 


(1) 


(2) 


(3) 


(4) 


(5) 


VLBA-PT 


17 


0.33 


0.76 


0.58 


VLBA-LA 


14 


0.47 


0.83 


0.69 


VLBA-FD 


13 


0.64 


0.90 


0.80 



Note. — Columns: (1) Site; (2) rms of the difference between consecutive pointings performed 
every 20 minutes during the observations; (3) loss of sensitivity due to the pointing errors: 
exp(— ln(2) error'^ /HWHM^) where HWHM is the half width at half maximum of the antenna beam; (4) 
loss assuming the telescopes are properly pointed at the beginning of each scan but drift to the pointing 
error values during a 20 minute scan; (5) loss given by a mixture of the two cases. 



Table 2: Description and Fit of Models for Sgr A* 



Model 


Vaj(™^) 


Vin(™^) 


PA (deg) 


x' 


X.' 




(7 


Comments 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


A 


0.34 ±0.14 


0.17 ±0.02 


22 ±20 


97 


1.19 







Best Fit Elliptical 


B 


0.18 ±0.02 






102 


1.25 


5 


1.5 


Best Fit Circular 


C 


0.175 


0.092 


80 


104 


1.28 


7 


2 


Scattering Alone 


D 


0.263 


0.175 


-10 


111 


1.36 


14 


3 


Lo et al. model scaled 


E 


0.44 


0.175 


-10 


247 


3.03 


150 


10 


Lo et al. model unsealed 


Note. - 


— Cohmms: (1) 


model name; 


(2) FWHM major axis; (3) FWHM 


minor 


axis; (4) position angle 



of major axis (for elliptical models); (5) of the model fit; (6) reduced obtained by dividing x^ by the 
number of baselines minus the rank of the matrix A^wA (see text); (7) x? difference from model A; (8) 
equivalent formal standard error with equivalent confidence at the Ax^ for the 3 degrees of freedom in an 
elliptical Gaussian model; (9) model description. 
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Fig. 1. — Comparison of three calibration methods for the VLBI array. Filled circles mark 
System Equivalent Flux Density (SEFD) measured using auto spectra on the hne calibrator 
VX Sgr. This method reflects the effects of pointing errors. The second calibration marked 
by crosses is the result of solving for station gains in the model fitting process described in 
the text. Similar degrees of scatter in the gains found by the two methods suggest they are 
both tracking similar variations. Pointing errors tabulated in Table |l| probably account for 
most of the observed gain variation. The third calibration used published gain curves and 
measured system temperatures to calculate the SEFD and is shown as solid lines. 

Fig. 2. — {u, v) coverage corresponding to all observations of Sgr A*. 

Fig. 3. — Calibrated visibility amplitudes for each baseline between Fort Davis (F), Kitt Peak 
(K), Owens Valley (O), Hat Creek (H), Pie Town (P) and Los Alamos (L). Calibration of 
each antenna was determined from a weighted least squares search through possible circular 
Gaussian models of brightness distribution. The dotted lines show amplitude curves for the 
best fit circular Gaussian model (FWHM=0.18 mas). Error bars are 3 a and the flux density 
scale for all baselines is the same. 

Fig. 4. — Calibrated visibility amplitudes as in Figure ^ plotted as a function of baseline 
length. The solid line shows the best fit circular Gaussian model and dotted lines show 
circular Gaussians that are ±1 a in size. Error bars are identical to those shown in Figure ^ 
and are left off here for clarity. The scatter in the data is reflective of the errors. 

Fig. 5. — vs the FWHM for a circular Gaussian model. 

Fig. 6. — Measured closure amplitudes for the Fort Davis, Kitt Peak, OVRO, Hat Creek 
quadrilateral along with closure amplitude curves for models A through E. The upper panel 
compares the closure amplitudes for models A through D. The lower panel shows model E 
and two variants: model E is the solid grey curve, the two dotted lines reflect 1 a errors in 
the N-S intrinsic structure at 43 GHz suggested by Lo et al. (1998). 

Fig. 7. — Size measurements of Sgr A* from VLBI at multiple wavelengths. Green squares are 
minor axis sizes and red dots are major axis sizes of best fit elliptical models at wavelengths 
less than 3.5 mm from Lo et al. (1998). Green and red solid lines are scattering relations 
given by Lo et al. (1998). The 3.5 mm upper limit to the size in the N-S direction from 
this work is shown as a horizontal green line. The N-S observed size predicted from the 
estimated 7 mm intrinsic structure (model D in Table 2) is the green filled circle, almost 
exactly at the upper size limit. Black triangles show sizes of circular Gaussian models fit to 
the 3.5 mm data described in this work, and to 1.4 mm VLBI data (Krichbaum et al. 1998). 
Hybrid ADAF model sizes from Ozel, Psaltis, and Narayan (2000) are added in quadrature 
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to both the minor and major axis scattering sizes and shown as the green and red dotted 
hnes respectively. The observed sizes predicted by the ADAF model exceed the 3.5 mm 
upper limits. 
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